First and second sound in a strongly interacting Fermi gas 
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Using a variational approach, we solve the equations of two-fluid hydrodynamics for a uniform and 
trapped Fermi gas at unitarity. In the uniform case, we find that the first and second sound modes 
are remarkably similar to those in superfluid Helium, a consequence of strong interactions. In the 
presence of harmonic trapping, first and second sound become degenerate at certain temperatures. 
At these points, second sound hybridizes with first sound and is strongly coupled with density 
fiuctuations, giving a promising way of observing second sound. We also discuss the possibility of 
exciting second sound by generating local heat perturbations. 
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I. INTRODUCTION 

First and second sound are a spectacular manifesta- 
tion of superfluidity in strongly interacting quantum liq- 
uids. They describe the coupled oscillations of the su- 
perfluid and normal fluid components at finite tempera- 
tures. In ^He, second sound is a temperature oscillation, 
in contrast to first sound, which is a density oscillation. 
That second sound should exist in superfluids was first 
pointed out by Tisza [Ij]. A complete theory based on 
the equations of two-fluid hydrodynamics was derived by 
Landau In conjunction with this theory, detailed 
measurements of the speed of second sound were crucial 
for developing a microscopic understanding of superfluid 
^He. Two-component Fermi gases close to unitarity 
provide a new type of strongly interacting superfluid. Be- 
ing a dilute gas, the unitary Fermi superfluid exhibits 
unique "universal" properties shared with neutron mat- 
ter and other exotic systems Q . An understanding of 
first and second sound in this system will be crucial to 
measuring quantities of interest such as transport coeffi- 
cients, the superfluid density, and the finite temperature 
equation of state. 

Many classic signatures of superfluidity have already 
been seen in ultracold gases [1, [j, 0] including quantized 
vortices, the absence of viscosity, and the Josephson ef- 
fect. Second sound has not yet been observed, however. 
Its detection requires that one be able to access the col- 
lisional hydrodynamic regime where the two-fluid equa- 
tions are valid. For Bose-condensed atomic gases, this is 
difficult because of the small value of the s-wave scatter- 
ing length. However, this requirement is easily satisfied 
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in superfluid Fermi gases close to unitarity where the 
s-wave scattering length is infinite. Due to the strong in- 
teractions characterizing the unitary regime, scattering 
between thermal excitations is sufficiently rapid to be in 
the hydrodynamic regime [1, [l^. At the same time, 
in strongly interacting atomic Fermi gases, second sound 
is weakly coupled to density perturbations. Modulating 
the frequency of the harmonic potential, for instance, will 
not appreciably excite second sound. Thus, even if exper- 
iments can reach the hydrodynamic regime, the ability to 
excite and detect second sound in trapped gases remains 
a significant challenge. 



In this article, we solve the Landau two-fluid equations 
of uniform and trapped Fermi superfluids at unitarity, 
presenting both numerical and analytical solutions. The 
nonuniformity of the equilibrium thermodynamic func- 
tions in a trapped gas makes a reliable solution of these 
equations very challenging. First solutions were obtained 
Ref. [ll[ using brute force methods with simplified ther- 
modynamic functions and in Ref. fl5| using a very simple 
ansatz for the velocity fields Vs{r,t) and v„(r,t). These 
papers did not, however, capture the basic features of the 
problem which represent the major achievement of the 
present paper: the surprising analogy between the uni- 
tary Fermi gas and superfluid Helium as a result of strong 
interactions, the strong dependence of second sound on 
the behaviour of the thermodynamic functions both at 
low T and close to Tc, and the peculiar bimodal struc- 
ture in the density response arising from the hybridiza- 
tion between first and second sound in a harmonic trap. 
In this work, we develop the variational approach to solve 
the two- fluid equations and illustrate these features. We 
also discuss the possibility of exciting second sound by 
generating local heat perturbations. 
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II. FIRST AND SECOND SOUND 

The dissipationless Landau two-fluid equations in a 
trap Voxt are given by UMM- 



dp 
dt 



ds 
dt 



V • (sv„ 



0, 



(1) 



at 



-VP-nVK=xt.(2) 



Here, j = psVs+p„v„ is the current density, ps and p„ are 
the superfluid and normal fluid densities for a gas with 
total mass density p = mn = Ps+Pn- P = —U + Ts-\- pm 
is the local pressure of the gas with energy density U and 
entropy density s, while p{n) is the local chemical poten- 
tial. Since we are interested only in the linear solutions 
of these hydrodynamic equations, we have omitted terms 
which are quadratic in velocity in ([2|). 

One can show [13, that the normal mode solutions 
of the hydrodynamic equations with frequency oj can be 
derived by minimizing the variational expression 



[5pf + 2p,[%)y-s+p,{^^)^{6-sf 



Jdr [psou2(r) -hp„o<(r)] 

(3) 

with respect to the displacement fields u„ and charac- 
terizing the density and entropy fluctuations according to 
6p = -V-[psQUs+p„oUn] and Ss = -u„- Vso-l-(so/po)V- 
[jOso(us — u„)]. These are related to the normal and su- 
perfluid velocity fields by = v^, u„ = v„. s = s/p 
is the entropy per unit mass. The effect of the trap- 
ping potential enters ^ through the position dependent 
equilibrium thermodynamic functions. These functions 
are calculated using the theory developed in Ref. p^ . 
which is an improved version [l6[ of the well known the- 
ory of Nozieres and Schmitt-Rink (NSR) ^ . We hope 
our work encourages more ab initio Monte-Carlo calcu- 
lations of thermodynamic functions [isj in order to im- 
prove the accuracy of the theoretical predictions for the 
propagation of second sound. 

Before discussing the general solutions of ([3]) , it is use- 
ful to solve this equation using ansatzes for a pure in- 
phase mode [Us(r,t) — u„(r, i) = u(^)(r, i)] and a pure 

(2) (2) 

out-of-phase mode [pso{r)ul ' {r,t) = -p„oU„ 0]> 
hereafter referred to as first and second sound. One can 
prove that these two modes satisfy an exact orthogonal- 
ity condition obeyed by the hydrodynamic equations (see 
Appendix |^ . These first and second sound modes cor- 
respond to pure density (6T{r, t) — 0) and temperature 
{Sp{r,t) = 0) oscillations, respectively fT2i|. They are 
the exact variational solutions of ^ when the coeffi- 
cient {dT/dp)s in ^ is set to zero. A nonzero value 
for this quantity leads to coupling between density and 
temperature oscillations and hence, between first and sec- 
ond sound. A convenient way of describing this cou- 
pling is in terms of the dimensionless Landau-Placzek 



(LP) ratio e 



where 



T{ds/dT)p and 



— T{ds/dT)p are the equilibrium specific heats per 
unit mass at constant pressure and density, respectively. 
When [dT/dp)s = 0, this means Cp = c„, or e = 0. 
We find that-similarly to liquid ^He-the solutions of the 
two-fluid equations for a trapped Fermi gas are well de- 
scribed by weakly coupled flrst and second sound modes, 
even though e can be order unity in Fermi gases. 

Inserting the ansatz Us(r, i) — u„(r, i) = vS^\r,t) 
for first sound into ([3]), taking the variation with re- 
spect to u'^"'^^ and making use of standard thermodynamic 
identities as well as the equilibrium conditions VPo = 
-noVl^ext and VTq = {dT/dp)sVpo + (ar/as)pVso = 
0, one obtains Euler's equation 



(1) 



-(V-u(i))VK, 



-V(u(i).Vl/cxt) 
(4) 



For a Fermi gas at unitarity, where P — 2U /i, it follows 
that po{dP/dp)s = 5Po/3, a result also satisfied by a 
noninteracting gas. Using this in one sees that it is 
identical to the equation for a coUisionally hydrodynamic 
Bose gas above Tc using ideal gas thermodynamics [H, 

(2) 

For second sound, we use the ansatz Us \r^t)pso{r) = 
—Un'{r,t)pno{r). Inserting this into the variational 
procedure gives the following equation for the superfluid 
displacement field: 



Po 



Po \ OsJ„ V PnO 



(5) 



Since Sp{r, T) — for this mode, one may also write 
this as a closed equation for the temperature fluctuations 
ST{r) = {dT / ds) p5s[r) , the entropy fluctuations 5s be- 

ing related to Us by the expression below ^ and the 
second sound ansatz. 



III. FIRST AND SECOND SOUND IN A 
UNIFORM FERMI GAS SUPERFLUID 

For a uniform superfluid (14xt — 0), the solutions of (|4|) 
and ([5]) are plane waves of wavevector q with dispersion 
LJi — Ciq and lu2 — C2q, where 



dP 



dp)-s' 



T 



si PsO 
Cv PnO 



(6) 



These flrst and second sound velocities are the standard 
expressions 0, [l^ used to describe superfluid ^He where 
e ^ 1 except in a narrow region around Tc. The mea- 
sured sound velocities, shown in FigJT]for superfluid ^He, 
as a function of temperature [20l. [2ll|7agree with (O. 

In Fig. [51 we plot the calculated values of ci and C2 
(shown by dashed lines) given by ^ for a uniform Fermi 
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FIG. 1: Experimental values of first and second sound speeds 
in superfluid *He. Calculating the thermodynamic functions 
in (O using the phonon-roton excitation spectrum measured 
by neutron scattering, one gets excellent agreement with this 
data (from Refs. 0|and [22]). 

superfluid gas at unitarity, using NSR thermodynam- 
ics [l^ m, [13 • In this figure, we also show the full 
solutions of the two-fluid equations, One immedi- 

ately notes that the full solutions are well described by 
the uncoupled first and second sound modes, ([6]). To un- 
derstand the effect of coupling (due to a finite value of 
{dT/dp)s in ([31)), we write the solution of the full two- 
fluid equations ([3|) as a linear combination of first (u*^^^) 

and second sound (us^„) modes. The resulting variational 
solutions of Q in a uniform superfluid are sound modes 
with dispersion [l^j 

^1,2- 2 ' ^ 

where 54 = Aq^4{pso/ Pno){dP/ds)l = 4^?^|e/(l+e) de- 
termines the coupling strength. (Note that {dP/ds)p = 
P[){dT I dp)s ) This shows that the coupling between 
first and second sound will be small as long as (c^ — 
c^Y /^c\c2 ^ e/(l -I- e), a condition met even if e is not 
small due to the fact that the speeds of first and second 
sound are never very close [23j . 

While the speed ci of first sound shown in Fig. [2| is 
weakly dependent on temperature, the speed C2 of sec- 
ond sound varies significantly between T = and Tc, 
where it vanishes (since p^ = 0). Comparing Figs. [1] 
and [H the similarity in the temperature dependence of 
C2{T) for a superfluid "^He and a uniform Fermi super- 
fluid gas at unitarity is quite striking. In both systems, 
Goldstone phonons are the dominant thermal excitations 
at low temperatures over a wide range of temperatures, 
which leads to the sharp increase in C2{T) as T ^ 0. 
Using ^ and the leading contributions to thermody- 
namic functions due to thermally excited phonons, we 
recover the celebrated results [1] ci = c and C2 — c/VS, 
where in our case c = y/S,/ Bt' p ~ 0.37vp is the T = 
value of the speed of sound [^. Our NSR theory re- 
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FIG. 2: Two-fluid sound speeds in a uniform Fermi gas at uni- 
tarity (solid lines). The uncoupled (e — 0) first and second 
sound speeds given by © are shown as dashed and dotted 
lines, respectively. Inset: Landau-Placzek ratio e as a func- 
tion of T/Tc. 



covers this T — result. If one instead used thermo- 
dynamics based on BCS Fermi excitations alone (with- 
out Goldstone phonons), one obtains the incorrect re- 
sult that C2 vanishes as T — > 0. The fact that phonons 
dominate thermodynamics up to such high temperatures 
(T ~ OATc) is unusual for a superfluid gas and results 
from the strong interactions exhibited by the unitary 
Fermi gas [IJl . These results only give the limiting low-T 
behavior of the two-fluid equations and arc not valid at 
very low T where hydrodynamics breaks down. At higher 
temperatures (T > 0.5Tc) thermodynamics is dominated 
by high energy rotons in superfluid ^He and by BCS 
quasiparticles (which have a large energy gap) in the uni- 
tary Fermi gas. The two-fluid modes discussed here are 
very different from those in dilute Bose-condensed gases 
where the weak interactions give rise to a strong cou- 
pling between density and entropy oscillations at all but 
the lowest temperatures [3, [l^l • 



IV. FIRST AND SECOND SOUND IN A 
TRAPPED FERMI GAS SUPERFLUID 

We now consider the case of a unitary Fermi gas con- 
fined by a harmonic trapping potential. For clarity, we 
only consider the breathing modes (Z = 0) for an isotropic 
potential Vcxt = ^mujQr^ . The extension to other types 
of modes (dipole, quadrupole, etc.) and to anisotropic 
potentials is straightforward in our approach. In Fig. [3l 
we show our results for the breathing mode frequencies 
of the full two-fluid equations using a variational polyno- 
mial ansatz (see Appendix IB)) for Us(r) and u„(r) in ([3]). 
Also shown are the frequencies of the lowest first [uji (n)] 
and second sound [w2('^)] modes, given by Q and ^ 
(n = 0, 1, ... is the number of radial nodes). First sound 
solutions are the "horizontal" branches while the dotted 
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lines are the second sound solutions. Figure [3^ a) shows 
that the frequencies of the uncoupled first and second 
sound modes are close to the full solutions of show- 
ing that, similar to the uniform case, the two-fluid modes 
of a trapped unitary Fermi gas are almost pure density 
and temperature oscillations. 

The two-fluid equations in ^ and ^ simplify to those 
of a pure superfluid {pno = 0) at T = and to a nor- 
mal fluid {pso — 0) above Tc. In both cases, there is 
a single fluid, whose displacement field is given by ([4]). 
In these two limits, (|H) gives analytic solutions for the 
first sound breathing mode frequencies. At T = 0, the 
hydrodynamic equation of motion for a trapped unitary 
Fermi superfluid [1^ predicts breathing modes with fre- 
quencies [1^, 113] uji (n) = luq ^y4:{n + l){n + 3)/3. Above 
Tc, our numerical results shown in Fig. [3l^a) for the first 
sound breathing mode frequencies quickly approach the 
analytic result for a coUisionally hydrodynamic classical 
gas 123], wi(n) = wo\/lOn/3 + 4. For the n = first 
sound mode, both expressions reduce to the same value 
wi(0) = 2luo- The n = first sound scaling solution 
cx r is actually an exact solution of the full two-fluid 
hydrodynamic equations at unitarity with this frequency 
at all temperatures pU. IT^ (moreover, it is an exact so- 
lution of the many-body Schrodinger equation at unitar- 
ity [11] for an isotropic trap). In contrast, the n = 1 
and higher first sound modes vary strongly with tem- 
perature. The measurement of the frequencies of these 
modes would be of great interest since they are sensitive 
to the temperature dependence of the equation of state 
of the Fermi gas. 

One can also derive a simple analytic results for the sec- 
ond sound frequencies uj2{n) as T ^ and also T Tc- 
For the low temperature region, our approach is anal- 
ogous to that used to obtain the uniform gas result 
C2 — c/VS. Assuming that phonons dominate thermody- 
namics, we solve ([5]) analytically within LDA. One finds 
fsee Appendix ICj) u!-y(n) = |(n-|-l)u;o. The numerical cal- 
culations of in Fig.[3|a) are difficult for T < OATc 
due to the smallness of the normal fluid density. However, 
we note that our low T results (T < 0.4Tc) extrapolate 
to our analytic results at T = 0. As noted in Sec. IIIH 
when T — > 0, the predictions of the two- fluid hydrody- 
namic equations are only useful in giving the limiting 
low-T behavior, the hydrodynamic regime for the ther- 
mal component no longer being ensured. 

As T ^ Tc, the behaviour of the second sound frequen- 
cies is very peculiar. Recall that the velocity of second 
sound in a uniform superfluid vanishes at Tc (see Fig. [2]) 
because p^o 0. The behaviour of second sound in a 
trap, on the other hand, depends sensitively on the spa- 
tial dependence of Pso(r) close to Tc. Within LDA, the 
size Rs = RTF\/i — T/Tc of the superfluid component 
vanishes as T Tc, giving rise to an increase in the min- 
imum wavevector, q ^ l/Rs- Using pso oc [1 — T/Tc]^^, 
we find that ^2 - C2q oc [1 - T/Tj^-i/^. For (3 < 1/2, 
LOi diverges. The p^ data we use |16| are consistent with 
the correct value /? ~ 1/3 [10] and consequently, we ex- 




FIG. 3: Left panel: Two-fluid modes of a unitary Fermi gas in 
an isotropic harmonic trap of frequency uio . The full solutions 
of the two-fluid equations are given by the solid lines. The 
lowest first sound modes (blue dashed lines) given by @ are 
shown along with the lowest second sound modes (red dotted 
lines) given by ([5]). The first sound mode frequencies ujiiji) 
are indistinguishable from the full solutions away from the 
hybridization points and have been shifted upwards slightly 
for clarity. The dashed lines at T = show our analytic 
results for the second sound frequency UJ2 (n) discussed in the 
text. Right panel: A blow-up of part of the left panel showing 
the hybridization between second sound and the n = 1 first 
sound mode. 

pect (and find) our calculated LDA values of ijJ2 diverge 
as r —*■ Tc. Our results emphasize that the low (T ~ 0) 
and high (T ~ Tc) temperature behaviour of the second 
sound mode frequencies depend sensitively on the ther- 
modynamic properties of the unitary Fermi gas. 

The rapid increase in the frequency of second sound as 
a function of temperature has the important consequence 
that, in contrast to the situation in a uniform Fermi gas 
(see Fig.[2|), the first and second sound mode frequencies 
cross each other. As anticipated in Ref. ^], this leads 
to a hybridization effect between first and second sound. 
Our results for the solution of the full two-fluid equa- 
tions ([2]) Jwhich includes the coupling between first and 
second sound left out of ^ and ([5])] clearly reveal the 
hybridization between the n — 1 and higher first sound 
modes and the second sound modes (being an exact so- 
lution of the full two-fluid hydrodynamic equations at all 
temperatures, the n = Q first sound mode does not couple 
to second sound). The hybridization between the n — 1 
first sound and the lowest four second sound modes is 
shown in Fig. EJb). These results are easily understood 
by looking for variational solutions of the full two-fluid 
equations ((31) as a superposition Ui.2 = AvS^^ + B\i)^)i 
of first and second sound. Varying ([3]) with respect to A 
and B, we obtain again the dispersion ([7]), where uji and 
ijj2 are the solutions of (|4]) and ([5]), while the coupling 
strength is now 

^2 ^ dv {dP/ds),V-{son^^)V-uW 
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Here, the integration is restricted to the superfluid region 

r<Rs and we define M(i)=i/dr[p,o(ui'^)^+p„o(ul;^)^]. 
The splitting between the two modes at hybridization, 
where wi = lj2, is given by — o)^ = (5^. As T — > Tc, 
this sphtting decreases (see Fig.[3][b)) since S vanishes at 
Tc where Rs 0. 



V. EXCITING AND DETECTING SECOND 
SOUND IN A TRAPPED GAS 

The results in Fig. [3l^b) show that the coupling be- 
tween second sound with first sound is strong in a Fermi 
superfiuid when their fi'equencies are equal. Hybridiza- 
tion means that second sound will couple to the density 
response at these crossing points. Adding a density per- 
turbation of the form 6V{r,t) = A/(r)e^"^* to the two- 
fluid equations, making the same ansatz as before for 
ui,2 (where now A and B depend on time), it is straight- 
forward to derive the response function [ij, for this 
perturbation, defined by / drf{r)Sp{r,t) = Ax(w)e~"^*, 
with imaginary part 

- = I' [ZiSi^' - ^l) + Z25{u^" - '^2)] • (9) 

Here, I = J drpou'^^ • V/(r)/-\/M(l) characterizes the 
strength of the density response to a perturbation with 
spatial dependence /(r). Equation ([9]) is vahd close 
to the hybridization points where the frequency of the 
first sound mode of interest crosses with the frequency 
of a second sound mode and we can ignore contribu- 
tions from other modes. The response is characterized 

~2_ 2 

by the two poles lui 2, with strengths Zi = rr-i — ^ and 
Z2 = 1 — Zi. At hybridization, where uji = L02, wc 
have Col 2 = ujf i and one finds Z\ = Z2 = 1/2. 
This result is extremely important from an experimen- 
tal standpoint. It shows, for instance, that switching 
on an appropriate static perturbation /(r) will result in 
a simultaneous excitation of the two modes, with the 
occurrence of a typical beating effect. The response is 
maximized by choosing /(r) such that V/(r) = u'^^^(r). 

Figured] shows the two fluid response, clearly showing 
the bimodal structure in the vicinity of hybridization be- 
tween the n = 1 flrst sound mode and second sound. For 
the n = 1 first sound mode, /(r) is well approximated 
by /(r) ~ — )^a{T)r^ where a = {r^) I {r^) is weakly 
temperature dependent. This potential can be straight- 
forwardly realized by modulating an optical dipole trap. 
Close to hybridization, this potential excites both first 
and second sound, as shown by Fig. d) Away from hy- 
bridization, it excites only the n — \ first sound mode. 
As discussed earlier, the measurement of this mode, char- 
acterized by the density fluctuation 8p ^ —V • [poV/], 
would be of great interest since it is sensitive to the equa- 
tion of state. 

Away from hybridization, the viscous damping of first 
sound should be small [3l[. Close to hybridization, on 




TIT, 



FIG. 4: Normalized density response Imx(ii') of a trapped 
Fermi gas at unitarity to a perturbation (5V^(r, t) = A/(r)e~'"' 
in the vicinity of the n = 1 first sound mode (see text). To 
simulate the effects of damping [3l| , the delta function peaks 
have been broadened by 0.026l)o- Inset: double peak structure 
at the hybridization point T ~ 0.64rc showing that first and 
second sound contribute equally (Z\ ~ Z-z) to the response. 
The frequency splitting Ao; ~ O.OBojq between peaks is large 
enough to be observed in experiments. 



the other hand, we expect that the damping will be en- 
hanced. The bulk viscosities and ^2 exactly vanish 
at unitarity [33| and the damping is determined by the 
shear viscosity 77, the bulk viscosity (^3, and the thermal 
conductivity k. When v,j = v„, the bulk viscosity C3 
and the thermal conductivity k (recall that 8T — ^ for 
this tj[pe of motion) do not contribute (see Sec. 140 in 
Ref. [231 and Ch. 19 in Ref. for discussions of hydro- 
dynamic damping). Thus, away from hybridization, only 
the small shear viscosity 77 contributes to hydrodynamic 
damping of first sound. In contrast, close to hybridiza- 
tion, where Vj 7^ v„, there will be enhanced damping 
arising from finite values of k and ■ If the damping is 
large enough, the double peak structure in Fig. [5] will be 
masked. However, the increased damping itself would be 
evidence of the excitation of second sound. 

A more direct way to excite second sound could be re- 
alized by generating a heat perturbation. For example, a 
laser tuned close to an optical resonance will heat the gas. 
The use of two intensity-modulated counter-propagating 
beams will produce a periodic heat perturbation, result- 
ing in an additional position and time dependent term in 
([T|) for the entropy density, leading to the excitation of 
second sound and generation of a large temperature fiuc- 
tuations when the frequency and the wavelength of the 
perturbation are in resonance with the second sound dis- 
persion. This procedure would be particularly efficient in 
an elongated trap where one could probe the dispersion 
relation ([S]) of a uniform superfluid [33|. 
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VI. CONCLUSIONS 

In this paper, we have presented reliable predictions 
for first and second sound in a Fermi gas at unitarity, in- 
cluding both numerical and analytic results. In contrast 
to the situation in weakly interacting Bose gases where 
the solutions of the two-fluid equations are, in general, 
a complex admixture of density and temperature oscil- 
lations (see, for instance, Refs. [H, [s^l), the two-fluid 
modes in a strongly interacting Fermi gas are weakly 
coupled density and temperature oscillations over a large 
range of temperatures. Nonetheless, we find that in the 
presence of trapping, the coupling can be significant in 
places where hybridization occurs between first and sec- 
ond sound. This effect provides a promising way to excite 
and detect second sound in trapped gases by measuring 
the response of the system to a density probe. Alter- 
nately, in highly anisotropic cigar-shaped traps, second 
sound waves could be excited directly by heating the gas 
with a space- and time-dependent perturbation or by a 
density perturbation (35| . 
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One immediately sees that as long as uji 7^ ^2, the exis- 
tence of a first sound solution (ui^^ = ul^') implies an or- 

(2) (2) 

thogonal second sound solution with Us Pso + u„ ' pno — 
0, and vice versa. 



APPENDIX B: POLYNOMIAL ANSATZ FOR 
THE DISPLACEMENT FIELDS 

Restricting ourselves to the breathing modes (/ — 0) in 
an isotropic trap, we use the following variational ansatz 
for the displacement fields: 

N-l 

u,(„) (r,t) = r a.(«)jr^'e-*"*, (Bl) 
3=0 

where N is the number of the variational parameters 
{cisj, o,n,j}- Inserting this ansatz into ([3]), the mode 
frequencies are obtained by minimizing the resulting ex- 
pression with respect to these N parameters. Fig. [3] is 
based on = 8. In actual calculations, we have found 
that it is better to work in terms of the variables Us — u„ 
and {psqUs + p„oU„)/po- 



APPENDIX A: ORTHOGONALITY CONDITION 

Suppose that {ui"'', ul"'-'} and {ui''\ul^^} are two so- 
lutions of ([3]). Inserting the first solution into the varia- 
tional equations that result from ([3]), left-multiplying by 
the second solution and integrating by parts, we find 

J dr psoui''^ ■ u(^) + iu] J dr p^ou^^^ ■ u^f ^ = 



dr 



-Po 



■ 1 




-po 




dT^ 


) Sp^ 

S 


dp, 



5(fc) 



(Al) 



Here, ujj is the frequency of the {ui"'-', u^"*-*} solution and 
are the corresponding density and entropy 
fiuct nations. 

We see that the right-hand side of Eq. (|A1[) is symmet- 
ric with respect to the exchange of the j and k solutions. 
Thus, we obtain the orthogonality relation 



drpsoui'^'^-ui^'^ + I drp„o<''^-uJi' 



,(fc).„(i) 



{u;^^u;l) = 0.(A2) 



APPENDIX C: LOW TEMPERATURE LIMIT OF 
THE SECOND SOUND FREQUENCIES IN A 
TRAPPED FERMI SUPERFLUID 



By combining the LDA formula Tpir) = Tf{Q)[1 - 
(r / Rtf)'^] for the Fermi temperature with standard ex- 
pressions for phonon thermodynamics [29| , we can reduce 
([5|) to an equation for temperature oscillations ST{r,t): 



{1 - f^)VlST + fdrST ^ 0. 



(CI) 



Here, we have used the dimensionless variables r = 
r/RxF and W2 = 2)L02IoJq. The result (jCip . which is valid 
in the low T two-fluid region, is similar in structure to the 
equation for the density oscillations in a Bose condensate 
at T = based on the Gross-Pitaevskii equation [3fi] . It 
can be solved exactly using analogous techniques. One 
finds that (5T(r, lo) ~ X)fc=o '^fc^^*^ are polynomial solu- 
tions with the eigenvalues (jJ2{n) = 2{n-\- 1). Thus, in the 
limit of low r, the breathing mode second sound tem- 
perature oscillations have the simple dispersion relation 
W2(n) = |(n + 1)ljo- 
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